Long-lived periodic revivals of coherence in an interacting Bose-Einstein condensate 
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We observe the coherence of an interacting two-component Bose-Einstein condensate (BEC) sur- 
viving for seconds in a trapped Ramsey interferometer. Mean-field driven collective oscillations of 
two components lead to periodic dephasing and rephasing of condensate wave functions with a slow 
decay of the interference fringe visibility. We apply spin echo synchronous with the self-rephasing of 
the condensate to reduce the influence of state-dependent atom losses, significantly enhancing the 
visibility up to 0.75 at the evolution time of 1.5 s. Mean-fleld theory consistently predicts higher 
visibility than experimentally observed values. We quantify the effects of classical and quantum 
noise and infer a coherence time of 2.8 s for a trapped condensate of 5.5 x lO** interacting atoms. 

PACS numbers: 03.75.Gg, 03.75.Dg, 67.85. Fg, 67.85.De 



Atom interferometry [J Q is a powerful method of 
precision measurements and a long phase accumulation 
time is desirable for improving sensitivity. Decoherence 
limits the timescale of an interferometric measurement 
and is of fundamental importance in quantum infor- 
mation processing. In this regard interparticle interac- 
tions can play a detrimental role in interferome- 
try with trapped Bose-Einstein condensates Non- 
linear interactions generate quantum phase diffusion Q 
and mean-field driven dephasing 0, [l^ which lead to 
the loss of interferometric contrast. Interaction-induced 
phase uncertainty has limited the coherence time in a 
multi-path interferometer to 20 ms [lH and in a double- 
well interferometer to 0.2 s Q. Monitoring the local 
spin coherence in the center of a two-component Bose- 
Einstein condensate (BEC) showed a promising decay 
time of 0.6 s [l2j . However, spatially nonuniform growth 
of the relative phase across the BEC leads to fast dephas- 
ing of the condensate order parameter and decay of the 
fringe visibility [lo| . 

Collisional dephasing can be reduced by tuning the s- 
wave scattering length to zero in the vicinity of a Fes- 
hbach resonance [l3|. Another way to minimise detri- 
mental interaction effects is to use noncondensed atoms 
with lower atomic density • A long coherence time of 
58 s was recently achieved in Ramsey interferometry with 
trapped cold atoms using rephasing via the identical spin 
rotation effect (ISRE) [iB] . 

In this Rapid Communication we demonstrate that the 
deleterious effect of atomic interactions on BEC coher- 
ence can be reversed via mean-field induced rephasing 
of the condensate wave functions. The periodic self- 
rephasing has a mechanism different from ISRE rephas- 
ing observed in noncondensed atoms and is due to 
induced collective oscillations of the condensate wave 
functions. The timely application of spin echo further 
enhances the visibility of Ramsey interference fringes 
for a trapped ®^Rb condensate and prolongs the cohcr- 
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FIG. 1. (Color) Simulations of collective oscillations of a 
two-component Bose-Einstein condensate: (a) axial density 
profiles of the components at different evolution times (solid 
lines) and the two-component ground state (dashed lines); (b) 
normalized local spin projection after Ramsey interferometry 
(iV = 5.5 X 10*, chemical potential 582 Hz, A/27r = -37 Hz, 
no losses). Lines of constant represent wavefronts of rela- 
tive phase. The mean-field dynamics act to periodically sep- 
arate the components and curve the phase wavefronts modu- 
lating the contrast of interference fringes. 



encc time to 2.8 s. The pseudospin-1/2 system consists 
of a two-component condensate trapped in two internal 
states |1) and |2), coupled via a two-photon microwave- 
radiofrequency (MW-RF) transition. Four factors af- 
fect coherence of the condensate: mean-field dephas- 
ing, asymmetric losses, quantum noise, and interaction 
with noncondensed atoms. With a single condensate 
produced in state |1), a 7r/2 pulse prepares a nonequi- 
librium coherent superposition of two states with mod- 
ified mean-field interactions and spatial modes distinct 
from those of the two-component ground state [Fig.[lja)]. 
Small differences in the s-wave scattering lengths an, 
ai2, and 022 cause the wave functions of the two com- 
ponents to oscillate out of phase around the ground- 
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state modes [16|. During the first half of the oscillation 
the relative phase grows inhomogeneously [l^l, as is evi- 
dent from the wavefront curvature of the normalized local 
spin density = (n2 — ni)/(ni + 71-2) [Fig. [TJb)], where 
ni = l^'ip is the density of state \i). In the second half 
of the oscillation, the wavefront curvature reverses and 
the relative phase becomes uniform along the condensate 
when the two components complete the cycle and the 
wave functions overlap again. The resulting self-revival 
of the interference contrast would be nearly complete if 
not for the presence of state-dependent loss of conden- 
sate atoms, which we compensate for by using a spin 



echo technique [17|, [1 



An *^^Rb BEC in the state IF = 1, tof -1) = |1) is 
produced on an atom chip [lOj in a cigar-shaped mag- 
netic trap with radial and axial frequencies 97.0(2) Hz 
and 11.69(3) Hz respectively. The magnetic field at 
the trap bottom is 3.228(5) G for which the first-order 
differential Zeeman shift between states |1) and |2) = 
|F = 2, mp = +\) is canceled [l^. This ensures that 
the relative phase is insensitive to magnetic fields and 
that atomic interactions dominate the dcphasing pro- 
cesses. Radiation for the two-photon coupling is provided 
by a half-wave dipole antenna (/mw ~ 6.8 GHz) and 
on-chip wires (/rf ~ 3.2 MHz). The fields are detuned 
by 1 MHz from the intermediate state |F = 2, = 0). 
The effective two-level system has a two-photon Rabi fre- 
quency r2i2/(27r) = 350 Hz and a two-photon detuning 
A/(27r) = -37 Hz. To detect both BEC components si- 
multaneously we use MW adiabatic passage to transfer 
99% of state |1) atoms to the \F = 2,tof = —1) state in 
free fall and apply a magnetic field gradient to spatially 
separate the components before taking a single absorp- 
tion image [lo| . In the time-of- flight images we observe 
no thermal fraction. 

The Ramsey sequence starts by applying a 7r/2 pulse 
of 0.7 ms duration, creating a coherent superposition of 
|1) and 1 2). Both components are confined in the mag- 
netic trap and after an evolution time T, an interro- 
gating 7r/2 pulse is applied. The phase lag Lp of this 
pulse is varied by a programmable MW source in order 
to obtain a Ramsey fringe in the phase domain. Imme- 
diately after the second pulse we release the atoms, mea- 
sure the populations iVi and N2 from an absorption im- 
age and evaluate the normalized atom number difference 
Pz = {N2 - Ni)/{Ni+ N2). In the mean-field formahsm 
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where N is the total atom number at T . To analyze the 
interference fringes we use the expression 

p,(iv,r,A,^) = v(iv,r)cos[Ar + 0(iv,r) + ^], (2) 

where 0(iV, T) is the interaction-induced relative phase, 
and V is the fringe visibility. Fringes in strongly de- 
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FIG. 2. (Color online) Ramsey interference fringes in the 
atom number [(a) and (c)] and phase [(b) and (d)] domains 
for evolution times of 20 ms [(a) and (b)] and 450 ms [(b) 
and (d)]. Run-to-run variations in the atom number A'^ cause 
correlated changes in data points for P^. This allows us to 
correct the phase of the data for the variations of A'' in phase 
domain (b) and (d). Solid lines represent expected fringes 
according to Eq. ([2}. 



pcnd on N as shown in Figs. [2ja) and[2jc). While ex- 
perimental control over A, T, and ip is straightforward, 
run-to-run fluctuations in the total number of condensed 
atoms lead to systematic shifts of the fringe. The ef- 
fect of such variations is corrected for by using the mea- 
sured number of atoms in each experimental cycle. We 
approximate (f>{N,T) = a{T)N'^/^T, where a{T) is de- 
termined from simulations of the Gross-Pitaevskii equa- 
tions (GPE) and slowly varies with T; q;(20 ms) = 0.79 
and q;(450 ms) = 0.90. This approximation for 4){N, T) 
agrees well with the fringes acquired by varying the atom 
number with constant phase tp [Figs.[2l^a) and [DJc)]. For 
each T we vary the microwave phase (p and obtain Ram- 
sey fringes in the phase domain ip -f- <f>{N,T), corrected 
for number fluctuations [Figs.lljb) and^Kd)], from which 
we evaluate the interference contrast V. 

In recent spin-squeezing experiments (20l - [2^ the atom 
number measurement was calibrated using the scaling 
between and the quantum projection noise. Wc cal- 
ibrate the atom number measurement using the simu- 
lated dependence of on the total number N at the 
evolution time of 340 ms (for 500 < iV < 10^). Wc nu- 
merically solve the three-dimensional coupled GPE equa- 
tions [IB,!!^ to obtain the wave functions ^'i and ^"2 and 
evaluate Pz using Eq. ([1]) . The best fit of the model yields 
an absolute calibration of N with a statistical uncertainty 
of 5%. 

In the temporal evolution of the Ramsey fringe visi- 
bility V{N,T) (initial total atom number 5.5(6) x 10"*, 
peak density 7.4 x 10^'^ cm~'^) we observe a decaying 
periodic structure with peaks and troughs [Fig. El[b)]. 
The initial decrease of visibility is due to the nonuni- 
form growth of the relative phase [l^ and the spatial 
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FIG. 3. (Color online) Schematic of interferometric sequence 
for (a) Ramsey and (c) spin echo. Corresponding plots of the 
visibility V [(b), (d)] where data (filled circles) are measured 
in the phase domain as a function of T. Curves correspond 
to simulations using coupled GPE theory (red dotted line), 
truncated Wigner model with classical noise (blue dashed), 
and with additional 10 s exponential decay (solid black). 



separation of the components [16[. By the end of the 
first collective oscillation the relative phase is uniform 
and the two components overlap again (Fig. [I}. This 
periodic evolution continues with a slow decay of fringe 
visibility (decay time of 1.3 s). Inelastic two-body col- 
lisions remove atoms from state |2) faster than from 
state |1) [2^ limiting the maximum fringe contrast to 
Vmax = 2y/W[N^/{Ni + N2). Thc couplcd GPE model 
[Fig. |3|b), dotted line] correctly predicts the period of 
the visibility revivals (Trev = 0.37 s), but overestimates 
their magnitudes. 

We apply a spin echo pulse [Fig. [3jc)] to reverse the 
relative phase evolution and compensate for asymmetric 
losses by inverting the populations. Maximum revivals 
of the visibility arc observed (V = 0.85 at T = 0.75 s 
and V = 0.75 at T = 1.5 s) when the sclf-rcphasing and 
spin echo work in unison. This occurs if the tt pulse 
is applied at T/2 = Ticv, 2 Trcv, ■ • when the relative 
phase is uniform along the condensate. The GPE simu- 
lations predict an almost complete revival of the visibility 
(V = 0.98 at T = 1.5 s), which is not seen in the experi- 
ments [Fig. [Sid)]. 

To characterize phase decoherence we studied the 
growth of the phase uncertainty for Ramsey and spin 
echo sequences (Fig. . For each time T we set the phase 
tf of thc second 7r/2 pulse to ensure (P^) = 0. In each re- 



alization we measure iV, record the deviation SPz.i from 
thc computed atom number dependence [Figs. [2ja) and 
Wis), solid line] and record the relative phase deviation 
5(j)i = arcsin((5Pj,i/V). We use M = 50 to 100 realiza- 
tions for each datum point in Fig. 2] and evaluate the 
Allan deviation a{(j)) for the duration of an experimental 
cycle (50 s). Error bars on each point in Fig. |4] represent 
thc statistical uncertainty of cr((?i) and arc calculated as 
a(0)/VM^. 

We evaluated contributions from three sources of clas- 
sical noise: MW frequency fiuctuations, uncertainty in 
the recorded atom number, and fluctuations of the split- 
ting (spin echo, recombination) process. To characterize 
the MW frequency instability we measured growth rates 
of thc phase uncertainty of 0.50(8) rad/s in the Ramsey 
and 0.125(20) rad/s in thc spin echo sequences (Fig. |4l 
dash-dotted lines) for 5.5 x 10^ uncondensed cold atoms 
trapped at 160 nK. Frequency fluctuations in our imag- 
ing laser contribute an uncertainty to the recorded atom 
number, increasing the spread of 5Pz,i [Fig. [SJa)]. The 
probe laser linewidth of 1.3 MHz (full width at half max- 
imum, FWHM) results in a 2.3% uncertainty in the mea- 
sured atom number when imaging on resonance. We find 
the contribution of the measured atom number uncer- 
tainty to the phase noise by calculating d(p/dN from the 
GPE simulations. 

We also considered the role of quantum noise in de- 
coherence of the condensate through quantum phase dif- 
fusion and nonlinear losses. We use the truncated 



Wigner approach [2J, |25[ and start from a Born-Markov 
master equation with additional terms for the domi- 
nant two-body losses The master equation is then 
transformed to a multimode Wigner representation, and 
thc resulting partial differential equation for thc Wigner 
quasiprobability function is truncated, leaving only first- 
and second-order derivatives. The equivalent stochas- 
tic equations are solved numerically, and the simulated 
quantum phase diffusion is plotted in Fig. |4] (dashed 
lines). In the experiment, the fluctuations of Pz af- 
ter splitting are not limited by quantum noise but are 
determined by shot-to-shot variations and drift of the 
two-photon Rabi frequency flu- Over several hours we 
typically measure a 2% deviation of (Pz) after thc first 
7r/2 pulse, which is above the standard quantum limit 
of 0.43%. After the inclusion of 2% splitting fluctua- 
tions into the simulations we obtain the phase diffusion 
induced by quantum and splitting noise. The black solid 
line (Fig. |4]) represents the combined contribution from 
the classical and quantum sources of noise and is very 
close to the spin echo data. 

The GPE simulations of the temporal evolution of vis- 
ibility (Fig. [3]) do not include irreversible decoherence 
effects and consistently show higher values than the ex- 
perimental data. If phase noise a{(f>) is the major source 
of decoherence the visibility would decrease by a coher- 
ence factor e~'^^'^^^^'^ [HI, [2^, which is equal to 0.92 for 
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FIG. 4. (Color online) Phase uncertainty growth in Ram- 
sey (a) and spin echo (b) sequences with the two-component 
BEC (black circles) and noncondensed cold atoms (red trian- 
gles) . The blue dashed lines show the effect of quantum noise 
simulated with the truncated Wigner method. The red dash- 
dotted lines indicate the contribution from MW frequency 
fluctuations. The black solid lines show the simulated com- 
bined effect of quantum and classical noise. 



In conclusion, we report that atomic interactions not 
only lead to quantum and spatial dephasing of a trapped 
BEC interferometer but also produce mean-field-drivcn 
rephasing through periodic collective oscillations. The 
synchronized application of spin echo significantly en- 
hances the self-rephasing, allowing us to record a visi- 
bility of 0.75 at an evolution time of 1.5 s. The quantum 
noise simulations consider the effects of asymmetric two- 
body losses and quantum phase diffusion on the BEC co- 
herence. Together with the evaluations of classical noise 
and single-particle coherence, the simulations allow us 
to estimate a condensate coherence time of 2.8 s, the 
longest coherence time for an interacting quantum de- 
generate gas. Future work will focus on the suppression 
of classical noise, reaching longer coherence times in the 
quantum noise dominated regime and exploring the effect 
of thermal atoms on condensate decoherence [J] . Our ob- 
servations show that long coherence times are possible in 
a trapped interacting BEC interferometer, which is an 
important result for precision measurements and sensor 
development based on quantum degenerate gases. 

We acknowledge helpful discussions with P. Hannaford 
and S. Whitlock. The project is supported by the ARC 
Centre of Excellence for Quantum- Atom Optics, and an 
ARC LIEF grant LE0668398. 



T = 0.75 s in the Ramsey sequence and for T = 1.5 s in 
the spin echo sequence; a value which cannot explain the 
detected decrease of visibility in both cases. In addition 
to the interaction-induced phase diffusion our quantum 
noise simulations take into account asymmetric two-body 
losses which are significant since 75% of condensed atoms 
are lost (predominantly through the state |2) decay) at 
the spin echo evolution time of 1.5 s. The modeling shows 
that at T = 1.5 s the mean value of the total pseudospin 
is (52)^/^ = 0.924(7V/2), an additional factor reducing 
the interference contrast. The combined effect of parti- 
cle losses and quantum and classical phase diffusion on 
the visibility in the spin echo and Ramsey sequences is 
shown in Fig. [3] (blue dashed line). In order to recon- 
cile the remaining deviation of experimental points from 
the simulation results we introduce single-particle deco- 
herence, possibly associated with the interaction of the 
condensate with uncondensed thermal atoms Q. A pre- 
vious study [l3| demonstrated that the condensate pres- 
ence shortens the coherence time of a thermal cloud, but 
the reverse effect of the uncondensed atoms has not been 
recorded. In our experiment the thermal fraction is small, 
as we do not observe its presence at the start or the end 
(1.5 s) of the interferometric sequence. We introduce an 
additional exponential decay with a time constant of 10 s 
to our simulations for the best fit with the spin echo data 
(black solid line in Fig. [3]) to infer an overall coherence 
time of 2.8 s. 
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